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Abstract 

We discuss three different ways to arrive at kaon condensation at Uc — 3no where 
no is nuclear matter density: (1) Fluctuating around the n = vacuum in chiral 
perturbation theory, (2) fluctuating around nvM near the chiral restoration density 
where the vector manifestation of hidden local symmetry is reached and (3) 
fluctuating around the Fermi liquid fixed point at ~ no- They all share one common 
theoretical basis, "hidden local symmetry." We argue that when the critical density 
Uc < Uy. is reached in a neutron star, the electrons turn into K~ mesons, which go 
into an S-wave Bose condensate. This reduces the pressure substantially and the 
neutron star goes into a black hole. Next we develop the argument that the collapse 
of a neutron star into a black hole takes place for a star of M ~ 1.5Mq. This means 
that Supernova 1987A had a black hole as result. We also show that two neutron 
stars in a binary have to be within 4% of each other in mass, for neutron stars 
sufficiently massive that they escape helium shell burning. For those that are so 
light that they do have helium shell burning, after a small correction for this they 
must be within 4% of each other in mass. Observations support the proximity in 
mass inside of a neutron star binary. The result of strangeness condensation is that 
there are ~ 5 times more low-mass black-hole, neutron-star binaries than double 
neutron-star binaries although the former are difficult to observe. 
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1 INTRODUCTION 



While the phase structure of hadronic matter at high temperature both below 
and above the chiral phase transition temperature is being mapped out by 
laboratory experiments with invaluable help from lattice QCD, the situation 
with high density is vastly different. There is little information about matter 
above the nuclear matter density Hq ~ 0.16 fm~^ from laboratory experiments 
and there is practically no guidance from QCD proper since lattice - the only 
non-perturbative QCD tool available - cannot handle high density and the 
reliable perturbative QCD approach can access density regime only at asymp- 
totic density where such novel phenomena like color-fiavor-locked supercon- 
ductivity can take place. To theorists' disappointments, though, this density 
regime may be totally irrelevant to nature. We have a complicated landscape 
of phases theoretically predicted near and above a chiral restoration n^^ which 
is as rich as the phase structure of water but they are all based on models, the 
reliability of which is uncertain given the paucity of experimental supports 
and lack of reliable theoretical control. Consequently the plethora of differ- 
ent scenarios for dense stellar systems such as neutron stars, black holes etc. 
available in the literature offer little guidance for understanding compact star 
physics. 

The aim of this paper is to exploit a systematic effective field theory framework 
from three different vantage points to describe with some confidence the phase 
structure of dense matter crucially relevant to the formation of stable compact 
stars and its implications on the population of neutron stars and light-mass 
black holes. These approaches, particularly the second and the third, rely on 
one common strategy that hadronic systems under extreme conditions can be 
accessed reliably by hidden local symmetry [1] that incorporates the scaling 
property of chiral symmetry [2]. Our reasoning will be backed by detailed 
analysis of astrophysical observations. 

This paper consists of two parts, one on an important phase change in hadronic 
physics and the other on collapse to black holes in astrophysics. Our aim is 
to bridge these seemingly disparate branches of physical phenomena. In the 
first part, we discuss recent developments on the most likely phase transition 
in hadronic matter, namely, kaon condensation, as density increases beyond 
the nuclear matter density to ~ 3 no- We will develop the thesis that this 
is the first - and perhaps the last - crucial phase change at high density 
that matters for the fate of compact stars, leaving wide-open the possibility of 
other forms of higher-density phases involving quark matter etc. In the second 
part, we give compelling arguments why astrophysical observations strongly 
support that neutron stars of mass greater than ~ I.'oMq (where Mq is the 
solar mass) cannot be stable, as a consequence of which any compact star 
more massive than the maximum stable mass must be in the form of a black 
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hole. These two developments will then be joined to arrive at the conclusion 
that kaon condensation at ~ Suq implies ~ 5 times more low-mass black-hole, 
neutron-star binaries than double neutron-star binaries. 



2 KAON CONDENSATION 



Although QCD cannot provide at present useful and quantitative information 
for hadronic interactions at high density relevant to the physics of compact 
stars, there are three vantage points at which we have available rehable effec- 
tive field theory tools to work with. The first is the matter-free T = n = 
vacuum about which fiuctuations can be described by effective chiral field 
theory. Here there is a wealth of experimental data to guide model building 
and extrapolating beyond the normal matter density. The second is the other 
extreme regime of high temperature and/or density at which chiral symmetry 
is supposed to be restored, namely, the "vector manifestation fixed point" in 
hidden local symmetry formulated by Harada and Yamawaki [1] characterized 
by the gauge coupling going to zero and the interactions between hadrons 
becoming weak. The third is in between the two limits, i.e., in the vicinity 
of nuclear matter density at which the Fermi-liquid fixed point is located. 
Here both theoretical and experimental information accumulated since many 
decades in nuclear physics offers both guidance and control near the nuclear 
matter density and indicates what needs to be taken into account beyond. The 
results of these three approaches giving the critical density ric ^ 3no have been 
reported before. What is new is the conceptual link between them provided by 
the recent development of chiral dynamics in hidden local symmetry theory 
recently supported by a holographic dual approach in string theory that gives 
an effective field theory of QCD in terms of (infinite towers of both) vector 
mesons and baryons and pions. 

We will approach the state of matter believed to be present in the interior of 
compact stars from these three vantage points and converge to the result that 
kaons must condense in the vicinity of Suq. This then leads us to propose that 
kaon condensation is the most likely phase transition to take place in hadronic 
matter in the density regime that is relevant to the fate of compact stars. Any 
other possible phases such as e.g., quark matter, color superconductivity etc. 
that can appear at higher densities will then be academic issues for stellar 
objects, although interesting from a purely theoretical point of view. As we 
shall outline, the metastable behavior of SN1987A which gave off neutrinos 
for ~ 12 seconds and then disappeared in all respects is symptomatic of a con- 
densate of negative charge kaons. In a number of papers we have outlined how 
evidence of a condensate which sends neutron stars into low-mass black holes, 
black holes not much more massive than neutron stars, is already present, and 
to us quite convincing, in the spectrum of masses in neutron-star binaries. 
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It seems highly hkely that the densities ~ 3no are reached in neutron stars. 
It is touch or go whether they were reached in the progenitor of SN1987A. 
However, they certainly are reached in the evolution of binary neutron stars 
in the case where the two progenitor giant masses are more than 4% different 
from each other in which case the first born neutron star finds itself in the red 
giant envelope of the companion giant as the latter evolves. This neutron star 
accretes matter from the evolving companion giant hypercritically until its 
mass is increased ~ 0.7 Mq, reaching densities enough for kaon condensation 
which sends it into a low-mass black hole. 

If, on the other hand, the two progenitor giants are within 4% of each other 
in mass, they will burn helium at the same time and avoid the red giant stage 
of the second giant, which will then end up as a neutron star. 

2.1 Kaon Condensation As Restoration Of Explicit Chiral Symmetry Break- 
ing 

As the first estimate of the critical density for kaon condensation, we approach 
dense hadronic matter from the T = n = vacuum. This can be efficiently 
done by resorting to Weinberg's "folk theorem" [3] using effective chiral field 
theory. The basic idea is to construct an EFT Lagrangian that takes into ac- 
count all relevant degrees of freedom at low energy and all pertinent symmetry 
constraints, in particular, chiral symmetry and do a systematic chiral pertur- 
bation theory calculation to as high an order as feasible. This procedure is by 
now fairly well established for pion-nucleon interactions as well as, to some 
degree, for finite nuclei and is being currently extended to nuclear matter [4]. 
The first approach made in this spirit is the work by Kaplan and Nelson [5] . In 
this subsection, we describe a toy description that while perhaps oversimpli- 
fied, nonetheless brings out the essential feature of the physics involved while 
ignoring what we consider to be inessential complications. The key idea here 
is that kaon condensation can be viewed as the restoration by baryon density 
of chiral symmetry explicitly broken by the strange quark mass. 

To formulate this idea in the most economical way, it was found to be sim- 
plest to take the V-spin projection of the Goldstone boson sphere onto the 
V-spin circle, composed of K" and a in the 5'[/(3) chiral Lagrangian [6]. The 
Hamiltonian for explicit K~ chiral symmetry breaking brings in the kaon con- 
tribution to the explicit chiral symmetry breaking in the nucleon T^kn and in 
the kaon mass m^. The Hamiltonian for explicit chiral symmetry breaking - 
which is of 0{p^) in the chiral power counting - is 

H^SB = ^kn{NN) cos + ^mlU'' sin^ 
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c^^KNiNN) [^-^-Pj +lrny^'0' (1) 

where the last expression is obtained for small fluctuations 9. Here we have 
taken mean fleld in the baryon sector. Dropping the term independent of 9, 
we find 



"K — ""K 



IWk ) 



The T^KN is 

{rriu + ms){N\uu + ss|iV) 
{rriu + md){N\uu + dd\N) 



-'KN 



Gt^N (3) 



where m„ ^ ^ are the current quark masses and Ut^n is the pion cr-term. Lattice 
calculations [7] give 

^ 389(14) MeV. (4) 

We can go one order higher in the "effective mass" of the kaon (2) by incor- 
porating what is called "range term" which turns out be important quantita- 
tively [8]. It amounts to changing T^kn to 



Sfjv = (1-0.37^^ E;,^ (5) 

\ J 

where ujk- is the (anti)kaon energy 



ujk- = + (6) 



In addition to the scalar attraction (which can be thought of as giving the 
kaon an effective mass) , there is the vector attraction given by the Tomozawa- 
Weinberg term [9] - which is the leading term, i.e., 0(p), in the chiral counting 



Vtw = -^ri ^ -60 MeV— (7) 

where tiq is nuclear matter density and we have taken the empirical value for 
the pion decay constant 
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Fig. 1. Projection onto the a-K plane. The angular variable represents fluctuation 
toward kaon mean field. 

~ 93 MeV. (8) 



In a neutron star the electron chemical potential /ig- , which in the progenitor 
of SN1987A was ~ 220 MeV [10] is the driving force towards kaon conden- 
sation [11]. Once the in-medium kaon mass is less than /ig then the 
degenerate electrons in the electron cloud in the neutron star can lower the 
energy by changing into a kaon condensate. The kaons are bosons, so they will 
immediately go into an S-state condensate, which lowers the pressure substan- 
tially from that of the degenerate electrons in the cloud in the neutron star. 
After this drop in pressure, the neutron star will go into a black hole in a 
light-crossing time. 

The above captures the calculation of kaon condensation as it stood for many 
years that indicated that kaons could condense at ric ^ 3no. The principal 
lesson we learn here is that baryon density plays a prominent role, namely, that 
the underlying mechanism for condensation is the restoration of the explicitly 
broken chiral symmetry by baryon density. As noted in [12], kaons condense 
because the s quark mass is neither too light nor too heavy. 

The most complete calculation of possible kaon condensation along the gen- 
eral line described above was carried out by Thorsson et al. [10] with three 
different "^kni expressed in other terms in their chiral perturbation theory. 
This calculation contained the medium dependence of the kaon mass through 
TiKN (see eq. (2)) but did not contain the medium dependence of /^r, the in- 
medium pion decay constant as required by HLS theory which we believe is 
indispensable as we will see below. 



In HLS treatment, this is connected to what is called "intrinsic background de- 
pendence" resulting from Wilsonian matching to QCD, the condensates of which 
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One may raise objections to the above simple treatment on the ground that 
several mechanisms that could be important are left out in the treatment. Most 
obvious are the roles of the A(1405) and hyperons, much discussed in the liter- 
ature, which when one does fluctuations from the matter-free vacuum, cannot 
be ignored. When looked at in terms of an RG flow, these degrees of freedom 
give "irrelevant" terms for the location of the critical density, that is, insignifi- 
cant for the precise value but they play a determinant role in triggering phase 
transitions. 

Since the role of the A(1405) in nuclei and dense matter has been greatly dis- 
puted in the past, we clarify what it as well as other related degrees of freedom 
- such as p-wave hyperon interactions - that we will refer to as "dangerously 
irrelevant operators" following condensed matter terminology, do in the con- 
densation phenomenon. The point we will drive at is that while they may play 
a crucial role in driving the system toward instability, once the direction in 
which the system moves is determined, the location of the critical density is 
highly insensitive to the strength of such interaction terms involved. 

We focus on the contribution of A(1405) to kaon-nuclear interactions. The 
A(1405)-nucleon interactions can be written as four-Fermi interactions which 
have canonical dimension 6. As such, it is "irrelevant" from the usual dimen- 
sional counting. To show what such a term does in meson condensation, we 
follow the argument given in [12] for generic meson condensation in strong 
interaction physics. It is a toy-model argument but it is applicable to the 
present case. The technique used is the one developed by Shankar [13] for 
Landau Fermi liquid theory. There is a distinct difference, however, between 
boson condensation and Fermi liquid in that while a "mass" term is relevant 
in the renormalization group sense, so figures importantly in the condensation 
process, the effective mass in the Fermi liquid is a fixed-point quantity, so it 
does not flow. This is because the Fermi surface in the latter is fixed by fiat. 

Let the meson field be denoted generically as $. We will focus on two terms, 
one "mass" term and the other interaction term, with the action written 
schematically (apart from the kinetic energy term) 



are background - temperature, density etc. - dependent. 

^ Similar arguments can be made for p-wave hyperon interactions. They give rise 
to "irrelevant" four- Fermi interactions. 



S — Sm + Si 



(9) 



where 




(10) 
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Fig. 2. Schematic diagram of RGE flow for the mass M. The dotted hne with arrow 
up shows the flow for D < |^ and the sohd line with arrow down for D > f"^- 

Si = j{dud?qY{ded'^kfh^*^^^^5\u,e,q,k) (11) 

where \1/ is the barvon field containing nucleons, hyperons and A(1405) etc. 
Of course both M Li| and h subsume many terms with varying powers in the 
suitable expansion (e.g., the chiral expansion in chiral perturbation theory). 
For our discussion we need not specify them here. 

The scaling law for the fields involved under scaling transformation deter- 
mined from the kinetic term (which is required to be invariant under the 
scaling) shows that the "mass" term is relevant but the interaction {h) term 
is irrelevant. Now the RGEs are found a la Shankar [13] by making the scale 
change A — > sA with < s < 1 where A is the cut-off scale [12], 



dM 

"df 

dh 



M-Dh, (12) 



- = -ft/2 (13) 

where t = — Ins and D = dn with d > is a constant depending on A/kp 
where kp is the Fermi momentum of the baryonic system. The solutions to 
the RGEs are 



^ Tilde represents that it contains other "mass-like" terms. 

^ For instance, h will contain, aside from the kaon mass, leading 0{p) term of 
~ ^K^KN^N, the next-to-leading order OijP') terms, ~ ^K^KN'^N etc. 
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M(t) = (Mo - lDho)e' + '^-Dhoe-"\ 
h{t)=hoe-'^\ 



(14) 
(15) 



As one decimates down with s — (t — > oo) integrating out higher lying 
modes, the boson-fermion interaction gets exponentially suppressed. We see 
that ii D < then the M flows toward +oo. On the other hand, for 



3 Mo 

the M is bound to turn negative, signalling instability against condensation. 
This can be seen in Figure 2. Now D — dn and d > 0, so that means that (16) 
can be satisfied by dialling the density. What matters is the "least irrelevant 
term" figuring in ho that drives flow toward condensation while more irrelevant 
terms - unless of course numerically enhanced - get suppressed. 

Once the system, flowing toward instability, reaches the critical point, then 
what matters will be the least "irrelevant" term that "eats up" the mass, with 
other irrelevant terms being suppressed. In fact a detailed analysis shows that 
the A(1405) that brings the initial attraction to direct the RG flow flgures 
negligibly in determining the location of the critical point [8] : it is found that 
change of the A(1405)-nucleon interaction strength by orders of magnitude 
affects the critical density by only a few %. 

Another objection that could be raised is that no account is made of the strong 
repulsion that can be generated in nuclear interactions in dense medium that 
would prevent kaon condensation at such a low density. In [14], Carlson ct 
al. argue that due account of short-range nuclear correlations could push the 
critical density above ~ 7no. This objection will be answered in what follows 
in terms of "hadronic freedom" in the vicinity of the critical density which 
prevents the repulsion from growing at high density. 



2.2 Strangeness Condensation By Expanding About the Fixed Point Of the 
Harada- Yamawaki Vector Manifestation 



As the second estimation of the critical density which we favor over the oth- 
ers, we calculate kaon condensation going top-down from the chiral restoration 
point. For this we adopt the hidden local symmetry approach to chiral dynam- 
ics developed by Harada and Yamawaki (HY) [1]. It is appealing to think of 
this theory as a truncated infinite-tower description of holographic dual QCD 
that arises from string theory (explained below) in which all vector mesons 
lying above the lowest ones, p, lo and 0, are integrated out. This theory is 
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characterized by that it has local gauge invariance, which allows a systematic 
chiral perturbation expansion with the vector mesons included and that the 
effective theory is matched to QCD a la Wilsonian RGE and hence "knows" 
about QCD. The crucial outcome of the HY's analysis is that the theory has 
what is called "vector manifestation (VM)" fixed point to which the theory 
flows subject to constraints imposed by QCD when the system is driven by 
either temperature [15], density [16] or the number of flavors [1] to the critical 
point at which chiral symmetry is restored. The flxed point is given by 

p*:=(^*,a*,/;) = (0,l,0) (17) 

where g is the hidden gauge couphng constant, a — F^/F^ is the ratio of 
"parametric decay constants" of the pion and the longitudinal component 
of the p meson F^. and is the physical pion decay constant. 

What is highly pertinent to us is that near the VM fixed point, fiuctuations 
of the hadrons figuring in the theory become well-defined, making processes 
taking place near the VM fixed point amenable to a simple calculation. It has 
been observed that even if a given process takes place away from the VM flxed 
point, so g ^ and /jr ^ 0, taking a ^ 1 can be a good starting point for 
describing fluctuations in certain processes. For example, a = 1 is an accurate 
approximation for describing the chiral doublers in heavy-light mesons such as 
D and D* [17], the 7r+-7r° mass difference [18], the nucleon EM form factors [19] 
etc. 

Assuming that kaons condense near, but below, the critical density for chiral 
restoration n^, we start from the VM flxed point (17) and make a renormalization- 
group flow analysis as density decreases from the VM flxed point as was done 
in [20]. In performing this calculation, we exploit two observations on: (1) the 
relevant degrees of freedom near the VM flxed point and (2) the weakness of 
interaction strength in what is referred to as "hadronic freedom" zone. 

We discuss these observations. 

The flrst and most important observation is that the relevant degrees of free- 
dom for kaon condensation are not those manifest in free space where the 
elementary kaon-nucleon interaction takes place, but more relevant are the 
degrees of freedom in the vicinity of the flxed point, since the condensation is 
in the neighborhood of this point at which chiral restoration takes place. 
Thus, for example, since {qq) is near zero, the HLS gauge coupling constant g is 
near zero, making the vector mesons nearly massless as pions are. This makes 
chiral perturbation calculation feasible and reliable with the vector mesons 
and pions treated on the same footing with their masses appearing as small 
parameters of 0{p^). Furthermore, the A(1405) and other excitations which 
are important near the matter-free vacuum, become irrelevant in the sense 
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defined above, so do not figure in the calculation of the critical density. 

We thus have vector mesons and pions as the explicit degrees of freedom. But 
then what about baryons? 

In holographic QCD which generalizes Harada-Yamawaki's HLS theory, there 
arc no fcrmion degrees of freedom. Baryons therefore must arise as solitons. 
Indeed in the holographic dual QCD model of Sakai and Sugimoto [21] which 
has correct features of chiral symmetry properties of QCD proper, baryons 
do arise naturally as instantons in five dimensions in which low-energy QCD 
is dualized or alternatively as skyrmions in four dimensions embedded in an 
infinite tower of vector mesons. The chiral properties of the resulting baryons, 
in particular, nucleons, are very well described in the approximations adopted 
by the model [22]. In order to exploit the VM fixed point, we need to integrate 
out all except the lowest members of the tower and treat hadonic matter in 
terms of skyrmions emerging from HY's HLS theory. Indeed hadronic matter 
at high density has been described in terms of skyrmion matter with an infinite 
winding number [23]. However there is an indication that skyrmions are not 
stable at some density n flash above Uq [24]. We will argue that this implies 
that the nucleon goes over to constituent (or quasi) quarks at n flash < ^x- 
We will first develop the argument in temperature and then apply it to the 
density case. 



2.3 Hadronic Freedom, the Flash Point And Why the Constituent Quark 
Model Is Relevant 

What happens when hadronic system is heated to near the critical temperature 
T^, i.e., from 120 to 175 MeV, was studied by Brown et al [25] using the STAR 
peripheral data [26] . For the peripheral experiments 120 MeV was the freezeout 
temperature; i.e., the hadrons leaving the system at that temperature did 
not interact further. The upper temperature T = 175 McV was taken to be 
the phase transition temperature at which the VM dictates that the HLS 
coupling g go to zero and the parameter a go to 1. Brown et al. [25] showed 
that the widths r(p — > 2tt) for p decay into two pions were zero at T = and 
remained small until just before freezeout at T = 120 MeV. This was because 
the F's went as (m*/mp)^ and the mass of the p was taken to increase from 
zero at to 90% of the on-shell F at 120 MeV. The latter value was called the 
"flash point" because the interactions of the p's at that temperature was that 
of (nearly) on-shell hadrons. The p-mesons in the experiment were recreated 
by following the pion tracks in the time projection chamber back to the vertex, 
identified as the p, for those which were emitted in p-decay. The 120 MeV was 
thus not only the freezeout temperature, but also the temperature at which 
the hadrons turned into constituent quarks. 
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We learn from the above observation that the region of temperature from 
~ 120 to thus dehneates a region in which interactions are neghgible, a 
region governed by the vanishing couphng constant g ~ {qq) 0. We can 
also deduce what happens to the a parameter which is defined by 



dominance of the pion EM form factor, a = 2 in matter-free space. In the 
large A^^^ limit in which the bare HLS Lagrangian is defined by Wilsonian 
matching, one finds a ~ 4/3 with the a = 1 point lying on the RG trajectory 
matched to QCD. The a = 2 point is not on that trajectory, implying that 
in HLS theory, a = 2 can only be an accident. Indeed as mentioned above, 
nature prefers a ~ 1 at least in the leading order of chiral perturbation theory. 
It has been shown by Harada and Sasaki [27] that in heat bath, the vector 
dominance given by a = 2 is strongly violated as a goes to 1 when temperature 
goes toward T^. This suggests that we associate a flash ~ 4/3 as the flash point 
at which T = T flash ~ 120 MeV. We interpret the region between a ~ 4/3 
at 120 MeV and a = 1 at to be a region in which the relevant fermionic 
degrees of freedom are constituent quarks, with the off-shell quarks flowing 
freely from higher to lower temperature. Nothing detectable goes on here in 
this region - it is just a region of space-time without anything except off- shell 
particles within it, hence called "hadronic freedom" region. 

In dense matter, the situation is much less clear, so we are going to make a cer- 
tain number of guesses and develop a scenario that parallels the temperature 
scenario. 

As mentioned, it is established that the vector manifestation fixed point (17) 
is reached when density is at the critical density n-^ at which (qq) = [16]. It 
has also been verified that rrip ^ g oc (qq) — in the vicinity of the critical 
density and (a— 1) oc {{qq)Y 0. However the question is: where is the density 
fiash point above which the system goes into the hadronic freedom region? In 
order to answer this question within the framework of HLS theory, we need the 
baryon degrees of freedom and at present, there is no systematic calculation 
in HLS theory that includes baryons. What we know now is described in [19]: 
Up to the nuclear matter density uq, the available phenomenology in nuclear 
processes indicates that the HLS coupling g stays more or less unchanged. 
Beyond no, we have no information from data but our basic assumption is that 
the fiash point n flash lies at the point from which g drops rapidly arriving at 
g = Q dX n = ric- From what we discussed above, we deduce that it is at n flash 

^ There is an indication from a skyrmion description of dense matter that a ~ 1 at 
the "pseudo-gap density" which is identified as the flash density [28]. This may be 
a better approximation than a ~ 4/3 adopted below. 
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that going top-down, the constituent quarks go on-shell forming baryons as 
the couphng constant becomes strong leaving the hadronic freedom region. 
Now if one extrapolates linearly from a = 1 at n^^, to a = 2 at n = 0, we arrive 
at a ~ 4/3 when the density reaches the density 



n flash -^riQ. (19) 

This is the density counterpart of Tfiash — 120 MeV. The hadronic freedom 
region is identified as the interval n flash "^^Hp -^n^. The putative repulsion 
one expects naively from the vantage point of the matter-free vacuum cannot 
survive within this region and would not provide the mechanism often invoked 
for stabilizing a neutron star at as large as ^ 2 Mq [29] If] 



We can understand this in terms of an effective chiral Lagrangian that incor- 
porates Brown-Rho scaling (or "intrinsic density dependence (IDD)" in HLS the- 
ory) [4] . Consider four-Fermi (baryon) interaction terms in the effective Lagrangian 
relevant to the vector-meson channels: 

CUmNf ^ -^{N^^Nf - ^{N^^rNf + ■■■. (20) 

i 

Here the asterisk denotes the IDD (intrinsic density dependence) in the parameters, 
i.e., Brown-Rho scaling, appearing in the EFT Lagrangian and the ellipsis stands 
for other terms that are allowed by chiral symmetry. Suppose now that the param- 
eters in the matrix elements of the four-Fermi interactions (20) are expanded in 
power of the density operator n = N'JqN. The resulting Lagrangian will contain 
six- and higher-Fermi fields and when mean field is taken, will give rise to density 
dependent parameters in the theory. This means that certain chirally symmetric 
n-body (with n > 2) interactions are subsumed in the IDD of the coefficients in 
(20) when truncated to the four- Fermi terms. Since the mean field using the La- 
grangian is equivalent to doing Landau Fermi liquid theory [4,30], certain important 
many-body effects are therefore encoded chirally symmetrically in the IDD of the 
parameters. 

In this formulation, the repulsion resides in the uj channel in (20). In Harada- 
Yamawaki HLS theory, the parametric mass of the p and uj (or the lowest member 
of the infinite tower in holographic QCD) drops as density increases whereas the 
vector coupling g remains more or less unchanged up to n ~ no and then decreases 
after the "flash density" n flash ~ 2no. This means the coefficient Ci,^ will increase 
as g^/m*^ with the falling mass up to, say, n ~ no or slightly above. This repulsion 
up to that density provides the mechanism for the saturation of nuclear matter. 
From n flashy the VM (vector manifestation) starts becoming operative, with the 
gauge coupling and the mass falling at the same rate so that the ratio 5*^/m* ^ will 
remain constant at higher density. The repulsion will cease accordingly. 
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2-4 Calculation Of the Critical Density ric For Kaon Condensation 



We now have all the pertinent degrees of freedom in the vicinity of n^. They 
are the constituent (or quasi) quarks, the (pseudo)Goldstone bosons vr and K 
and the vector mesons p and u with their masses comparable to the pion mass. 
To describe kaon condensation, we need to take into consideration interactions 
between the kaon and the quasiquarks. The Lagrangian we have is HLS La- 
grangian that is matched to QCD at a density n ^ at a matching scale 
A\.^. Now we need to calculate RGE flow in theory with the gauge coupling g 
and the parameter a as described above in the density region n flash ^ '^x- 
The only other quantity that we have to consider is the parametric pion decay 
constant F*, since kaons will couple to the quarks with the constant l/F*. 
We need this constant at the matching scale A^. For this we observe that F* 
does not scale appreciably in the matter-free space as found by Harada and 
Yamawaki [1]. We shall give an argument why we expect this to be similar in 
medium. For this, it will suffice to to know what it is near the nuclear mat- 
ter density. This information comes from the analysis of deeply bound pionic 
atoms [3ljZ] from which we learned that the in medium parametric F^ must 
be decreased 



~ 20% at n = no- We have no systematic calculation of F* beyond the nuclear 
matter density. This is because F* runs with the intrinsic density dependence 
dictated by the Wilsonian matching of the HLS correlators to those of QCD, 
which is unknown since lattice techniques cannot handle density effects. How- 
ever one can argue that it will stay more or less unchanged up to the critical 
density. In both temperature and density, when one goes beyond the flash 
point, the vector meson mass rriy and the gauge coupling constant g* scale 
in the same way, i.e. proportionally to (qq)*- Hence the ratio g* /my l/F* 
- with a* ~ 1 - will stay more or less unchanged. This feature was exploited 
in the high-temperature case in [33]. As for the density relevant to kaon con- 
densation, we need not go all the way to the chiral restoration point but one 
can make a rough estimate of the ratio F^(n^)/-F^(0) in a HLS model with 
quasiquarks [16], 



It thus seems justified to take the scaling (21) operative from the fiash point 
to the kaon condensation point. It was shown in [20] that the Walecka vector 

^ A detailed review of the pionic atom data and the evidence for F* = 0.8F^ is 
given by E. Friedman and A. Gal [32]. 



F^ ^ F:{no) ^ rAno) ^ 0.8F^ 



(21) 




(22) 
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mean field used so much in nuclear physics already had this increase in it, so 
that it is part of the usual phenomenology of nuclear physics. 

Since the gauge coupling constant g is near zero in the hadronic freedom 
region, quantum loop corrections can be ignored. It suffices therefore to do 
the mean field calculation (i.e., tree contributions). For the neutron fraction 
Xn and proton fraction Xp in a compact star, the potential felt by the kaon 
will then be 



We have included p as well as uj exchange so as to make a typical neutron star 
charge up at density n^. Note that (23) is not in linear-density approximation. 
In this picture, this is the only important contribution to the effective mass 
of the kaon. Other terms are suppressed by the hadronic freedom effect. Now 
let us see how the effective kaon mass behaves near the chiral phase transition 
density riy, ~ Auq- For this we take a* — 1 and F*/Ft^ ~ 0.8. Then we find 



\9v / '^'v ] fixed point 



/my] zero density [o*-?"^^] fixed point 0-8^ 



zero density ^ J_ ^ ^ (24) 



For Xn ~ 0.9 corresponding to the neutron-star matter, we have from (23) 



VK-{n = 4no) ^ -516 MeV. (25) 

This suggests that the kaon must have a vanishing mass at the chiral transition 
density. 

Let us move down to the flash point ~ 3no. At this density, taking into account 
the flow of a, we find 



rriK^^n^)- ^ rriK-/^ ^ 165 MeV. (26) 

This is somewhat smaller but not too far from the electron chemical potential 
expected at that density, ~ 220 MeV [10]. Considering the roughness of the 
estimate, it seems reasonable to identify the flash point as the kaon conden- 
sation critical density n^. We note that were it not for the presence of the 
electron Fermi sea, the flow of the star matter would end up at the VM fixed 
point at riy^. However, the decay of electrons into kaons stops the hidden local 
symmetry (HLS) flow. 
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2. 5 Kaon Condensation From the Fermi-Liquid Fixed Point 



The third way to arrive at the kaon condensation critical density ~ 3no is 
to fluctuate from the density at which nuclear matter is saturated, namely, 
no ~ 0.16 fm^"^. The basic idea has been presented in previous review articles 
by the authors, the most recent of which is given in [30]. As discussed there, 
going above the nuclear matter density consists of first writing an effective 
Lagrangian that describes nuclear matter as the Fermi-liquid fixed point a la 
Shankar when the mean field approximation is taken. Given such a Lagrangian 
with the parameters determined with the "intrinsic density dependence" as 
dictated by the Wilsonian matching to QCD prescribed in HLS theory, one can 
make a fluctuation around the Fermi liquid flxed point to go to higher density. 
The effective Lagrangian which is anchored on chiral symmetry is equivalent to 
Walecka's linear mean fleld Lagrangian with the intrinsic density dependence 
(i.e. Brown-Rho scaling) suitably taken into account. 

To proceed to kaon condensation, we assume as above that the relevant fermionic 
degrees of freedom are the constituent quarks. We consider the constituent 
quark picture applicable at nuclear matter density and above and that the 
nucleon is made of loosely bound three quasiquarks and the kaon of loosely 
bound strange quark and chiral anti-quasiquark. Now the mean fleld nuclear 
potential has two parts, the vector potential Vn and the scalar potential S'at, 
which can be viewed as due to exchange of an uj meson and a scalar meson 
(denoted here s). 



9 

V/v = 9 n, 

*2 

SN^-^,n. (27) 

Here we are introducing a flctitious scalar in the line of Walecka model with 
scaling mass and coupling constant but one can actually avoid the scalar 
degree of freedom by writing the effective Lagrangian in terms of a four-Fermi 
interaction (see, e.g., [30]). Now from the phenomenology of Walecka model, 
we know that 



SN{no) - VN{no) ^ -600 MeV. (28) 

Taking into account of the fact that the kaon contains only one non-strange 
quark compared with three in the nucleon, we expect the kaon to feel the 
attractive potential 



Sk- + Vk- ^ -200 MeV. (29) 
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Now applying the same reasoning to the compact star matter with a 90% 
neutron fraction, we find at n ~ 2no that 



ujK-{2nQ) leOMeV. (30) 

Since the electron chemical potential expected at this density is 173 MeV [10], 
we see that kaon condensation is expected to take place already at 2no in this 
scenario. Thus this picture predicts ric < 3no consistent with the two other 
estimates. 



2. 6 Delayed Collapse Of Neutron Star 



It is not generally realized that the metastabihty of SN1987A, i.e. the fact 
that it emitted neutrinos for about 12 seconds before disappearing (and was 

"never heard from again"), was probably aided substantially by the presence 
of strange particles in the core [34]. For this, the fact that the particles were 
strange was not the important point; the strange particles with strange quarks, 
be they the K~ or the hyperon E", have negative charges and additional pro- 
tons would be needed to neutralize their charge. The additional protons would 
interact strongly attractively with the neutron main component in the neu- 
tron star and increase the binding energy of the star. As the star collapses to a 
neutron star, once the densities are ~ 10^^ g cm~^ the neutrinos are trapped. 
That is, they may try to random walk their way out, in the frame of the in- 
faUing matter, but once the density is ~ 10^^ g cm~^ the mean free path for 
neutrino scattering off the nucleons is only a few centimeters, so the neutrino 
drift (random walk) velocity outwards becomes smaller than the velocity of 
infalling matter. The neutrinos, above the trapping density, equilibrate with 
the neutrons and protons so that 

/^n = /^p + /^e - l^u, (31) 

and the neutrinos fill a Fermi sea of Fermi energy ~ 100 MeV. Thus, the 
neutrinos contribute to the pressure. 

It has been considered likely that the neutron star has hyperons in it. 
Although the mass of the E~ hyperon is 1197 MeV, much heavier than the 
neutron 940 MeV, the neutrons will be degenerate. Since the admixture of S~ 
will have large negative charge, its presence would bring down the electron 
chemical potential /Xe considerably [36]. The E~ at rest can replace both an 
electron and a neutron. It begins to come in at a density such that 

/^S- = /^n + /^e- (32) 
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This is generally satisfied at a density n ~ 2no; i.e., at about twice nuclear 
matter density [35]. (The A's come in earlier, at a lower density, but the A is 
charge neutral and does not have the interesting effects that the S~ has.) As 
the neutrinos leave, the electrons recombine with the protons to make neutrons 
and E~'s, and the EOS stiffens as the number of protons decreases. Thus, the 
star will not collapse immediately into a black hole. Kaon condensation will 
set in later as the star cools down further. 

As has been discussed extensively in the literature, the appearance of hyperons 
in dense matter has been invoked to argue against kaon condensation setting 
in at a density ^ 3rio because the driving force, i.e., the electron chemical 
potential, would get weakened [36]. We believe this argument is not valid. We 
have made clear that nucleons and hyperons are no longer the relevant degrees 
of freedom when the flash point, which we estimate to be at ~ 3no is reached. 
The RG flow is run completely by the Weinberg-Tomozawa vector interaction 
here, the hyperons having been integrated out, as irrelevant. Our fluctuation 
about the fixed point rather than about the perturbative vacuum has given a 
scaling which increases the vector interaction by (a*)~^, increasing from 0.5 in 
free space to 1 at the fixed point. This explains much of the increased attraction 
and, therefore, softening of the EOS in our nonperturbative calculation. 



3 A SCENARIO FOR A LARGE NUMBER OF LOW-MASS BLACK 
HOLES IN THE GALAXY 



This was the title of a paper written by G.E. Brown and H.A. Bethe in W.K. 
Kellogg Radiation Laboratory [37]. Since that time, much improvements on 
the calculations have been made, the most important one of which being the 
expanding about the VM fixed point of hidden local symmetry. And further- 
more more observational data have been accumulated. 

In this second half of the paper, we would like to present a chain of astrophys- 
ical observations to link the kaon condensation in hadronic matter developed 
in the first half to the population of neutron stars and black holes. The most 
pertinent and daring prediction on this issue had been made by Brown [38], 
who had found the accepted scenario for binary neutron star evolution to be 
completely wrong. As we shall outline, the standard scenario accepted up to 
then assumed that the accretion onto the first-born neutron star, while the 
companion giant evolved and the neutron star and giant went into common 
envelope, was limited by the Eddington rate and, therefore, negligible. It was 
noted in [37], after the suggestion by Chevalier [39], that the rate of accretion 
necessitated going over to hypercritical accretion and that that was sufficient 
to send the first-born neutron star into a (low-mass) black hole. It was then 
suggested that the only way to avoid this was for the two progenitor giants 
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to burn helium at the same time, which allowed the neutron star to avoid the 
red giant stage of the companion giant. This point is further developed below. 

The Brown scenario had the prediction that the two neutron stars in a binary, 
which must have burned helium at the same time, should be within 4% of each 
other in mass in order to do this. As we shall see, this prediction is fulfilled 
by the most massive neutron star binaries. 

In the case of the less massive double neutron stars, i.e., J0737-3039 and 
J1756— 2251, the pulsars acquire an additional 0.1 — O.2M0 during the helium 
red giant stage of the companion, which does not take place in the more 
massive binaries, but aside from this, the pulsar and companion must be within 
4% of each other in mass. The companion in the less massive double neutron 
stars, i.e., J0737-3039 and J1756-2251 have to burn hotter than those in the 
more massive ones during the helium burning in order to reach the same central 
temperature of the more massive ones because the surface to volume ratio is 
larger in the less massive stars. This means that more energy is lost from the 
surface. The result of this is that the progenitors of these two less massive 
stars go through a helium red giant stage as well as (earlier) a hydrogen red 
giant stage. During the He red giant stage the pulsar accretes 0.1 — O.2M0 
of helium. Therefore, the necessary limit of at most 4% difference between 
progenitor masses of the pulsar and companion can be somewhat larger if the 
companion of the pulsar gives an additional 0.1 to O.2M0 of helium during the 
helium shell burning. Note that the helium burning of pulsar and companion 
progenitors must overlap in time. The helium red giant stage is put into the 
evolution of J0737-3039 by Willems and Kalogera [40]. 

We will try to make our arguments for the general physics community. The 
general reader would benefit from references to a few review papers, where 
neutron star structure [41,42,43], supernova explosions [44,45,46] and binary 
evolution [47,48] are discussed. 



3.1 Hypercritical Accretion Onto First-Bom Neutron Star During Common 
Envelope Evolution 

In the study of neutron stars, a key tool used in astrophysics is the acceptance 
of Eddington accretion as the maximum possible rate. Matter falling on a 
neutron star is bound to the neutron star gravitationally by ~ 20% of its 
rest mass because of the great mass of the neutron star. This binding energy 
is used to produce radiation, say in the form of photons, which are emitted 
from the surface of the neutron star. These photons scatter off the incoming 
matter; which we take to be hydrogen, and deposit their momentum in it, by 
Thomson scattering. At the Eddington limit 
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MEdd = i?6 1-5 X 10-«MNsyr-' 



(33) 



where Rq is the radius of the neutron star in units of 10^ cm (the latter being 
a typical neutron star radius) and Mns is the mass of the neutron star in 
units of the mass of our sun Mq, the outward pressure from the scattering of 
photons by the hydrogen is sufficient to counteract the inward gravitational 
force. This, then, gives a maximum rate of accretion in the neighborhood of 
rates ~ MEdd- ^Edd is the rate of mass accretion which just holds off the 
accreting protons, 

MEdd = 2 X 10^^ ergs/sec. (34) 

(We shall see later on that accretion can begin again at M ~ lO^MEdd, where 
it is called "hypercritical accretion.") 

Now the processes that we are going to discuss such as common envelope 
evolution take a time ~ years. Given the 10~^ in Eq. (33), accretion at the 
order of the Eddington rate is unimportant and can be neglected. 

Let us begin our evolution of binary neutron stars starting from two giants, 
with zero age main sequence (ZAMS) masses 16 and 12 Mq. Each star burns 
hydrogen for ~ 90% of its lifetime, but the lifetime of the 16Mq star is shorter 
than that of the 12Mq. When the former has finished core burning of hydrogen, 
it expands in red giant stage, and the outer part of the red giant crosses the 
Roche Lobe, the line of equipotential between the two stars. The more massive 
16Mq star evolves first into red giant transferring ~ 8Mq to 12Mq star by 
Roche Lobe overfiow leaving the He core and 2OM0 star as in Fig. 3. 

In binary evolution usually a fiat q {q = M2/M1) distribution in the IMF 
(initial mass function) is assumed. This has the correlation between primary 
and secondary star in it that the former is more massive than the latter, 
which must be present becaiisc it evolves first. Although offhand this looks 
like a quite minimal correlation, it produces factor ~ 2 difference in results. 

We want to zero in on the common envelope evolution of the first born neutron 
star with the helium core of the 20Mq star. The common envelope evolution 
takes only ~ years, so thinking of Eddington limit for mass accretion of 1.5 x 
10^*^1/0 yr^^ the accretion is quite negligible and has generally been neglected 
in the literature to date. However it is unreasonable to assume the Eddington 
rate to hold forever, as the rate of accretion continuously increases. Ultimately 
the density of matter falling onto the neutron star must be sufficiently high 
that the photons originating from the neutron star by random walk in the 
framework of the infalling matter actually move out at a lower velocity than 
the infalling matter has in falling onto the neutron star, so that the photons 
are swept back onto the neutron star by the adiabatic infiow. Good discussions 
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WMc. star UM^. star 




Roche Lobe overflow 




He core 20M(7) star 




Collapse of He star into a neutron star 




20 Mq star evolves into red giant, enveloping the neutron star 

Fig. 3. Conventional model for double neutron star binary evolution. The neutron 
star goes into common envelope evolution, dropping towards the star, the increased 
gravitational binding energy being used to expel the hydrogen envelope. The He 
star explodes into a neutron star and if the two neutron stars remain bound, this is 
a binary pulsar. 
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of this are given by [49,50,51,52]. 



With the increase in the rate of accretion, which we scale with Eddington 



M 



m — 



Edd 



(35) 



in the region of rates m ~ 1, the cutting off of increase in accretion above that 

value sets in. The accretion from m ~ 1 to 10^ may be complicated, the effects 
from the Eddington limit causing uneven or sporadic flow, but for m > 10^ 
where hypercritical accretion sets in, the accretion again becomes simple with 
uniform flow. 



The nature may not be so simple because the matter is accreted at a low 
temperature, which means that it cannot get rid of its energy rapidly by 
neutrino emission, which increases with a high power of temperature T, but 
is low for the accretion temperatures. 

In the scenario of Bethe and Brown [56], neutrinos were trapped in the collapse 
of large stars, so the trapping of photons in order to give hypercritical accretion 
seemed self-evident. Indeed, one can work out that as the accretion increases 
to the Eddington limit, the photon trapping radius moves to the radius of the 
neutron star it^NS ~ 10^ cm. There is, then, a range of M's, 



MEdd <M< lO^MEdd (36) 

in which the outcome of accretion is not clear, some matter being accepted 
by the neutron star and some matter being lost in space. This dees not really 
concern us because for rates ~ lO'^MEdd and higher, the accreted matter piles 
up in an accretion shock. As the amount of matter increases, the density at the 
base of the accretion shock increases and ultimately, when the temperature 
reaches T ~ 1 MeV neutrino pairs carry off the energy as rapidly as it is 
accreted, holding the temperature to ~ 1 MeV. An analytical calculation of 
this was carried out by Brown and Weingartner [50]. 

Given the two giant progenitors going through the evolution shown in Fig. 3, 
and the helium stars, the lowest mass in the problem was the neutron star 
mass Mns, and the equations could be solved algebraically if Mns was set 
equal to zero [56]. The result was that the flrst-born neutron star (see Fig. 3) 
accreted IMq in the common envelope evolution during the red giant phase 
of the remaining giant. Belczynski et al. [53] (denoted as BKB) removed the 
approximation made by Bethe and Brown and obtained the result that the 
accretion was ~ O.SMq for the lower mass neutron stars and ~ O.QM© for 
the higher mass ones. In either case it was sufficient to send the primary into 



23 



a black hole (see Case D2 of BKB). Zeldovich et al. [54] had found already 
in 1972 that the accreting matter piles up, pushing its way up through the 
accretion disk. At this stage we go over to the two-dimensional hydro cal- 
culations of Armitage and Livio [55]. They showed that an accretion disk 
reformed inside of the accretion shock, allowing matter to accrete onto the 
neutron star. They, however, suggested that jets might drive the hypercritical 
accreting matter off, saving the neutron star from going into a black hole. 
We find, however, in agreement with observations, that the first born neutron 
star (the pulsar) is within 4% difference from the mass of the companion in 
the binary, in agreement with our scenario that hypercritical accretion can be 
avoided if the two giant progenitors are within 4% of each other in mass. No 
sign of any other amount of hypercritical accretion, which could simply add 
to the mass of the neutron star is seen. If, indeed, the accretion were stopped 
at some intermediate stage by jets, one would expect a distribution in pulsar 
masses compared with companion masses that differed from the 4%. Since no 
such distribution is seen, our assumption that the pulsar has gone into a black 
hole (if the giant progenitors differ by more than 4% in mass) is validated. 

Bethe and Brown [56] found that the accretion was ~ 10^ to 10^ times Ed- 
dington accretion, and assumed that with this vast amount of matter, one was 
back to the classical accretion of Bondi [57] which took no note of the above 
complication, but did guide the matter correctly through the sonic point; i.e., 
enforced the classical requirements for accretion. The accretion rate could not 
be determined directly in the Bethe and Brown [56] calculation, but it was 
compared with the fully three-dimensional numerical calculation of Terman et 
al. [58] 



3.2 Formation Of Double Neutron Star Binaries 

We do see double neutron star binaries, so the first-born neutron star cannot 
always go into a black hole. If the two giant progenitors are so close in mass 
that they burn helium at the same time, then they would go into a helium 
common envelope evolution, the one star sending He across the Roche Lobe 
into the other and the other one sending its expanding hehum to the first 
star [38]. Helium burning takes up 10% of the star lifetime, most of the time 
being spent in the main sequence hydrogen burning. To go from lifetimes to 
masses one must divide by about 2.5, so the two giant progenitors must be 
within 4% of each other in mass. 

This was a definite prediction modulo one assumption: That helium is not 
accreted during the helium common envelope evolution. The assumption was 
justified however by Braun and Langer [59] who showed that the helium burn- 
ing time is too short for appreciable mass to be accepted by either helium star. 
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Table 1 

Compilation of NS-NS binaries [60] . 



Object Mass (Mq) Companion Mass {Mq) References 



J1518+49 


-'-•00_Q44 




[61] 


B1534+12 


1 qqqo+0.0010 

1.000Z_Q QQJ^Q 


i.040Z_o 0010 


[61] 


B1913+16 


i.44U8_0.ooo3 


1 0070+0.0003 

l.OO/O_0.0003 


[61] 


B2127+11C 


1 Q4Q+0.040 


1 q«Q+0040 
i.ouo_o 040 


[61] 


J0737-3039A 


1 qc!7+0.005 
-'-•oo/_ooo5 


1.250l[!:[j{]^ (J0737-3039B) 


[62] 


J1756-2251 




1 1O+0.03 
-•-•-^0-0.02 


[63] 



Thus, the clear prediction is that two neutron stars in a binary must be within 
4% of each other in mass. 

There is, however, one proviso. That is that the lower mass neutron stars in the 
double pulsars J0737— 3039 and J1756— 2251 go through not only the hydrogen 
burning red giant but also a helium burning red giant. The temperature for 
helium burning in low mass stars must be greater than in higher mass stars in 
order to ensure a sufficiently high central temperature because of the energy 
loss through the surface. During the helium burning red giant, ~ 0.1 to O.2M0 
can be deposited on the first born neutron star by the helium star companion, 
and the first born neutron star should be that much more massive than the 
other, in addition to the possible ~ 4% difference in mass because they must 
burn helium at the same time. 

In Table 1 wc show the compilation of masses of double neutron star binaries. 
In J1518+49 the uncertainties are too large to say anything, except that the 
largest error +0.45Mq is in the direction of equahzing the masses. In the 
next three binaries the masses are very close. In fact, B2127+11 is in the 
globular cluster and often considered to have resulted from the exchange of 
neutron stars formed in different binaries. However, they do satisfy our nearly 
equal mass scenario. In the evolutionary calculations of the double pulsar 
J0737— 3039A,B the hehum red giant phase has been put in [64,65]. The 0.1 — 
O.2M0 accretion during the He red giant phase was calculated in hypercritical 
accretion by Lee et al. [66]. Using a flat distribution for the Initial Mass 
Function, Lee et al. calculated in hypercritical accretion what the masses of 
primary and secondary neutron stars would be in the scenario shown in Fig. 4 
in which case the primary (first-born) neutron star is assumed not to go into 
a black hole, but the neutron star of mass it ends up with after accretion is 
assumed to be stable. One sees that the primary star, the pulsar, has high 
~ 84% probability of ending up with mass between 1.8 and 2.3M0. This is 
not seen at all in Table 1, so we assume that these must all have gone into 
black holes. 
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Fig. 4. Masses of primary compact stars with and without accretion during H red 
giant stage of secondary star [66] . The filled circles show the primary (more massive) 
pulsar masses in double neutron star binaries. The corresponding secondary ZAMS 
masses were obtained from the secondary (less massive) pulsar masses using eq. (8) 
of Lee et al. [66]. Note that the 84% corresponds to Mcompact Star,primary > l.SM©. 
There is uncertainty in the final primary compact star masses due to the extra mass 
accretion, ~ O.IM© — 0.2Mq, during He giant stage. This may increase the primary 
compact star masses for both 'with-' and 'without-accrction'. Note that with our 
maximum neutron star mass of I.SMq, all primary compact stars with accretion 
would go into low-mass black holes. 

Note that only the primary pulsar masses (the first born pulsar in the dou- 
ble pulsar) would be affected by the accretion, the companion masses lying 
on the lower thin band in Fig. 4. In the figure, we have added the primary 
pulsar masses of 5 double neutron star binaries in filled circles. Wc left out 
J1518+49 because of the large uncertainty in its masses. We see that the a 
priori probability for the primary pulsar masses to lie on the lower thin band 
within 0.2Mq uncertainty due to the possible accretion during the He giant 
phase, given no other possibility than to remain neutron stars, is 



Obviously there is negligible probability assuming random neutron star masses 
that the two neutron stars in each of 5 binaries are all within 4% of each other 
in mass. 



P = (0.16)^ ~ 1 X 10-^ 



(37) 
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Tabic 2 

Neutron star masses in neutron-star, white dwarf binaries in which the central value 
of the neutron star mass is > 1.5Mq [67]. 



Name Pulsar Mass (M( 



J1713+0747 l-54;o;?o^ 

B1855+09 ^■57toil 

J0751+1807 2.10lg|g 

J1804-2718 < 1.70 

J1012+5307 l-68lg:i 

J0621+1002 l-70ig:29 

J0437-4715 l-58l[]:i^ 

J2019+2425 < 1.51 



For our later argument, we note here that the Hulse-Taylor pulsar B1913-I-16 
has mass 1.44O8±O.OOO3M0. It is the most massive - and the most accurately 
measured - of the neutron stars in binaries. 



3.3 Neutron Star - White Dwarf Binaries 



In Table 2, we list the neutron-star, white-dwarf binaries in which the central 
value of the neutron star mass is > 1.5Mq. We should point out from the 
neutron star masses in Table 2, given the errors noted, that there is no neutron 
star mass other than in J0751-I-1807, which is definitely above 1.6Mq. The 
interval from l.SM© to 2.IM0 of 0.6Mq is about half of the total spread, from 
I.IM0 to 2.IM0 in measured neutron star masses, so it would appear strange 
for only one of them to lie in the upper half of the possible interval. 

On the other hand, sufficient of the central values of the neutron stars in 
binaries with white dwarf lie above I.5M0 that it might be prudent to quote 
the maximum neutron star mass as the Bethe and Brown [68] value of I.56M0 
deduced from the cut in the Fe production in 1987A, rather than the earlier 
I.5M0 [37]. 

The conclusion of Lee et al. [66] is that the maximum neutron star mass must 
be less than I.8M0, otherwise the spectrum of double neutron star binaries 
would be completely different from that shown in Table 1. 

The largest measured mass is in J0751-I-1807, a neutron star in a binary with 
white dwarf, with mass 2.ltoioMQ [69]. However at 95% confidence level, the 
mass is 2.11q;5M0 so it could be as low as I.6M0 at this level. Clearly it 
could be I.8M0 with not unreasonable probability. As mentioned below, the 
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measured value has recently been reduced to 1.26Io;i2M0. 

A number of calculations were made to estimate the mass of the compact 
object in SN1987A. They were all more or less in the region of masses arrived 
at by Bethe and Brown [68] from the ~ O.OTSM© of Ni which was ejected in 
SN1987A. This gave an estimate of the separation in the density, between the 
mass that fell back onto the compact object and that which was ejected. The 
deduced compact object mass was < 1.56Mq. Calculations by Thielemann 
et al. [70] estimated the mass of the compact object formed in SN1987A to 
be within the range of gravitational masses 1.30 — 1.50 Mq. The higher val- 
ues follow from accretion during the formation of the delayed shock. The 
delayed scenario is now believed to be the correct one. The 1.5Mq pertaining 
to SN1987A is supported by the work of Timmes et al. [71]. 

Tauris and Savonije [72,73] evolved neutron star white-dwarf binaries, mostly 
by conservative mass transfer, so the evolution should be reliable. Most of 
their neutron star masses were in the neighborhood of 2.1Mq whereas all of 
the other 8 accurately measured neutron star masses in binaries with white 
dwarfs [60] are close to 1.5Mq. There must be some reason that limits them 
but we do not know what it is at the moment. 

We are well aware that neutron stars in binaries in some globular clusters 
have been reported with masses substantially above our I.56M0 limit. For 

example, the relativistic periastron advance for the two eccentric systems in 
the globular cluster Terzan 5, 1 and J, indicates that at least one of these 
pulsars has a mass of 1.68Mq at 95% confidence [74], not so different from the 
I.6M0 of the Nice et al. J0751+1807 I.GMq at 95% confidence. 

The situation in globular clusters is very different from Galactic. With very 
low mctallicity, their stars resemble more Population 11 stars than Galactic 
ones. This should not make a difference however in our calculation of kaon 
condensation and its consequences since we find the condensation depends 
only on the baryon number density. 

Our statement on kaon condensation is admittedly a very strong one. But it 
is a falsifiable prediction and should be straightforward to prove us wrong: find 
a well measured double neutron star binary in which the two neutron stars are 
more than 4 % different from each other ( modulo some small additional shift 
by He red giant) in mass. 

3.4 Are There Two Branches of Neutron Stars ? 

If the 2.IM0 neutron star of Nice et al. [69] were firmly established, then there 
would have to be a loop-hole or an alternative to the Brown-Bethe scenario 
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for the formation of neutron stars and light-mass black holes. Indeed Haensel 
et al. [75] have proposed an interesting way to reconcile the Brown and Bethe 
M^g*^ = 1.5Mq with the Nice et al. [69] 2.1Mq neutron star, two branches 
of neutron stars — reconciling a 2Mq pulsar and SN1987A. Two scenarios of 
neutron star formation were considered. Here we interpret the two scenarios in 
terms of two possibilities in our language. In the first scenario the formation 
of the strangeness core is rapid; this is the case of the double neutron star 
binaries in which the hypercritical accretion from the evolving companion 
onto the first-born neutron star, the pulsar, is at a rate of ~ lO^MEdd, so that 
the total accretion time is ~ 1 year. In the second scenario the neutron star is 
already present when the main sequence progenitor of the white dwarf begins 
to evolve. The neutron star accretes matter (and is spun up) by accretion from 
the companion red giant in stable mass transfer. The rate of mass transfer is 
^ MEdd, so it is ~ 10^ times slower than in the first case. Models of quark 
deconfinement (producing quark matter) are found in which the condition 

The Haensel et al. reconciliation is ingenious and deserves to be paid at- 
tention. However in view of the recent development on Nice et al's 2.1Mq 
neutron star discussed in the next subsection, reconciliation of that type does 
not seem called for. We would nonetheless suggest that the reasoning for the 
two tracks of Haensel et al is not valid for the case at hand. Of the two pos- 
sible phases, quark matter and kaon condensation, the former phase would 
be favored at high temperature, as is found in collapses like 1987A, and the 
latter in a cool scenario like the ~ 10* yr accretion from an evolving main 
sequence star onto the neutron star to make a neutron-star, white-dwarf bi- 
nary. In 1987A the collapsing star heats up, to a temperature ~ 25 MeV in 
the center, because of Joule heating. The latter is accomplished by the ~ 100 
MeV highly degenerate neutrinos leaving their degeneracy energy in the cen- 
ter of the forming neutron star as they leave with energies determined by the 
surface temperature. On the other hand, the temperature from accretion onto 
neutron stars in the white-dwarf, neutron-star binary evolution is much less 
than 1 MeV, so the thermal neutrino emission is almost negligible. Thus, the 
collapse into strangeness condensation would be favored by the more massive 
of the two stars in a double-track scenario and the transition into a black hole 
by strangeness condensation in the collapse of a large star would favor quark 
matter. In other words, the two tracks of Haensel et al would be filled in the 
different way, if indeed they were. 



3.5 New Measurement Of the Neutron Star Mass in J0751+1807 

A new value for the mass of J0751-I-1807 announced by D. Nice in the McGill 
meeting on "40 Years of Pulsars," August 12-17, 2007, http: / / www.ns2007.org^ 
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is now 1.26loj2J^0 to replace the old value (2.1 ± 0.2)Mq. This removes the 
necessity for two branches of neutron stars, with which we had troubles noted 
just above and reinforces the argument given in [76] that the 2.1 Mq must 
be wrong, for were it correct, the spectrum of neutron-star masses would be 
completely different from what it is. Thus, our statement is that all ends with 
baryon number density ~ Suq, there is only "nothingness," the space-time 
volume of the black hole, beyond. Assuming that we are correct, this leads 
us to the strong statement that color superconductivity, color-flavor locking 
and all of the other "exotic" phases that have been suggested for high density 
matter lie in the "black-hole nothingness," being inaccessible by the observers. 

3.6 Neutron Stars In Globular Clusters 

We have no doubt that neutron stars with masses above our limit of Mats (max) = 
I.56M0 will continue to be reported for globular clusters. No real evolutions 
of these have been carried out. From the measurement of the rate of advance 
of the angle of periastron of the binaries combined with matching white dwarfs 
with the data base, a Bayesian-type analysis gives the possibilities of different 
white-dwarf masses. Together with the mass function, probabilities of neutron 
star masses are then obtained. 

Binaries with measurable quantities, i.e., those with rapid advances in the rate 
of periastron advances, are clearly favored and the probability of measurability 
should be built into the Bayesian analysis. 

Why we were so worried about J0751+1807 was that the Shapiro effect was 
involved in the analysis. With an accurate Shapiro effect and mass function, 
the two masses in the binary can be measured accurately. The Shapiro effect 
is a relativistic effect resulting in the time retardation of radiation from the 
companion passing by the neutron star. Thus, it can only be measured when 
the angle of incidence is ^ 80°, so that the neutron star and companion are 
nearly in the same plane. The measurement of this shift gives a relation be- 
tween that of the two masses additional to that of the mass function, so with 
the two relations one can obtain both masses. 

Stars in globular clusters have generally low metallicity, more like Population 
II and III stars. They tend to be much more massive than Galactic stars. Our 
assertion that kaon condensation with subsequent collapse into black holes 
depends only on the baryon number density will be very interesting to test in 
the new environment of globular clusters. Following the present preliminary 

^ They will surely generate - as the 2.1Mq mass did - a large number of theoretical 
publications that will purport to invalidate the onset of kaon condensation at ~ 3no 
and its crucial role in the fate of compact stars. 
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"Schadenfreude" of observers trying to prove that theorists are wrong we look 
forward to real tests of our assertion in environments with low metallicity. 



3. 7 Convective Carbon Burning And Compact Star Masses 



Although the star in the center of SN1987A disappeared, we do know that 
the original star Sanduleak-69°202 was in star catalogs and did have a mass 
of ~ 18Mq. For many years there was substantial uncertainty in the evolution 
of giants in the region of masses in which they would begin, with increasing 
ZAMS mass, to go into black holes. The most important nuclear process in 
this evolution was the 



^2C(«,7)^^0 (38) 

reaction. This reaction determines where the convective burning of carbon 
stops. In many years of work, the Stuttgart group [77] and more recently 
Buchmann and Barnes [78] determined this to be 



f 162 ± 39 keV barns (Stuttgart) 

[ SQtlo + ^^-18 keV barns (Buchmann & Barnes) 

The upper superscript 300 means that its value is given for a reaction energy 
of 300 keV. The convective ^^C burning carries off a lot of entropy through 
neutrino pairs; it is the first reaction in the stellar burning of stars to do 
so. Once the ^^C is sufficiently used up so that this reaction no longer takes 
place, the increase in entropy must come from the increase in size of the Fe 
core, which has entropy ~ 3/4 per nucleon (a bit less than the ~ unity of the 
original Bethe et al. paper [79]). We show in Fig. 5 that the convective ^^C 
burning ends at I8M0 [80]. Note that the Fe core mass increases above 1.5Mq 
just at ~ I8M0. (Brown et al. point out that the Fe core mass should be a 
good indicator of the compact object mass because the decrease in mass from 
the general relativistic effect is compensated for by fallback in formation of 
the compact object.) Thus, with the best available nuclear physics to date, the 
supernova progenitor of SN1987A collapsed into a compact object of ~ 1.5Mq, 
not far from the mass constructed by Bethe and Brown from the 0.075Mq of 
Ni formed in the explosion of < 1.56Mq. 



^ The standard reply of astronomers to the statement that if a neutron star were 
present in the region where SN1987A exploded, then one should have seen something 
has been: "The absence of evidence is not evidence of absence." We do not believe 
this to be true because one would have seen something had a neutron star been 
there. It could not have just disappeared, leaving no evidence. 
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Fig. 5. Comparison of the iron core masses at the time of iron core implosion for a 
finely spaced grid of stellar masses [81]. The circular black dots were calculated with 
the Woosley & Weaver code [82], whereas the crosses employ the vastly improved 
Langanke, Martinez-Pinedo rates [83] for electron capture and beta decay. If the as- 
sembled core mass is greater than Mpc = I.8M0, where Mpc is the proto-compact 
star mass as defined by Brown & Bethe [37], there is no stability and no bounce; 
the core collapses into a high mass BH. M^s = 1.5Mq denotes the maximum mass 
of NS [37]. The mass of the heaviest known well-measured pulsar, PSR B1913-I-16, 
is also indicated with dashed horizontal line [84]. 

Note that the low-mass black holes will only be formed from ZAMS ~ 18 — 
2OM0 stars; i.e., black holes of 1.5 — 2.5Mq. In fact, the black hole from 
1987A is our only example of such a black hole. The next lowest mass one is 
the ~ 5Mq black hole in the transient source GRO J1655 which comes from 
an ~ 3OM0 giant [85]. We would expect black holes from giants with ZAMS 
masses 20 — 30Mq to appear in gamma ray bursts, one of ~ 3Mq from ZAMS 
mass ~ 2OM0 in GRB060218/SN2006aj [86]. 
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Fig. 6. Distribution of pR^ vs included mass in the presupernova model of a star of 
mass I8M0, according to Woosley [87] 

3.8 Accretion Onto Neutron Star After Supernova Explosion 

The shock wave powering the supernova explosion moves outwards through the 
helium core, essentially as a Sedov self-similar shock of constant velocity, the 
pressure being essentially inversely proportional to the volume. However pR^, 
where p is the density and R is the radial distance from the center increase 
greatly in the hydrogen envelope (see Fig. 6). The shock has to slow down and 
the material behind it likewise [87]. 

The shock wave in the supernova after the first four seconds will be slowed 
down as it has to pick up more and more material in the mantle and envelope. 
The outgoing material in the shock wave thus has to be decelerated. Within 
this material, there is a density discontinuity where the He mantle adjoins the 
H envelope. In the deceleration, the lighter H gas has to push the heavier He, 
the classical situation for RT instability. 

The deceleration takes place because the product pr^ increases going out. The 
main increase, from 7 x 10^^ to 1.2 x 10'^'^, takes place continuously in the H 
envelope, from — 7 — IIM© (sec Fig. 6). In the Sedov theory, a shock wave 
continues at constant speed if pr^ =const because its velocity is 



and p is inversely proportional to the enclosed volume. But if pr^ increases, the 
shock has to slow down, and the material behind it likewise. This deceleration 
gradually sharpens into an ingoing shock. 



U ~ (p/p)'/' 



(40) 
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Woosley's calculation [88] shows that when the outgoing shock enters the 
hydrogen envelope, the deceleration of the shock sharpens into a reverse shock 
and this shock has been used to deposit ~ O.IM© back on the compact object 
[89] . Fig. 3 of Shigeyama et al. [90] shows the matter being deposited back on 
the neutron star. 

When material sent radially inward by the reverse shock reaches the supernova 
center, the compact object is no longer there in case the neutron star receives a 
kick velocity. The material proceeds to move outwards from the original center. 
Eventually all of the matter will be moving radially outward, through with 
different speeds depending on the distance from the center. Thus Chevalier 
[91] suggests that the compact object will end up traveling at the same speed 
as its kick velocity. Therefore, we may, to a good approximation, consider 
the compact object to be at rest with respect to its ambient matter and 
use the Bondi [57] spherical accretion theory to determine M, the rate of 
mass accretion onto the compact object. Brown and Weingartner [50] find 
that the photons from the initial accretion are trapped by the high density 
infalling matter, but that in ~ 0.6 yr one would have accretion at basically 
the Eddington limit. 



3.9 Observability Premium For Unseen Black Hole - Neutron Star Binaries 



Our prediction on the implication of kaon condensation on neutron stars in a 
binary has important ramifications on the number of light-mass black holes. 
We claim that there should be ^ 5 times more low-mass black-hole, neutron- 
star binaries than double neutron-star binaries. 

Now what is the status of the observation? 

The observed pulsars have all been "recycled". That is, they have accreted 
mass from their companion. Empirically, this accretion of mass [92] brings 
down the magnetic field of the pulsar. The time that a pulsar can be seen is 
inversely proportional to its magnetic field. If the magnetic field is brought 
down a factor of ~ 100 from the recycling, the pulsar will be observable 
for 100 times longer. It is convenient to introduce the idea of "observability 
premium" [93] 



n = ^ (41, 

where B is the magnetic field of the pulsar. 11 ~ 1 for fresh, unrecycled pul- 
sars, which would appear in neutron star, black hole binaries, because the black 
hole cannot transfer mass ("recycle") the neutron star. For PSR 1913-1-16, the 
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Hulsc-Taylor binary, B ~ 10^*^ G for the pulsar, so it will be observable for 
~ 100 times longer than the black hole, neutron star binary. In Table 1 we 
have 6 double neutron star binaries. Taking 100 to be a typical difference be- 
tween the observable times for recycled and unrecycled pulsars, we see that we 
should see ~ 0.6 neutron-star, black-hole binaries. In fact, the factor is prob- 
ably somewhat less than this because the observer would be somewhat biased 
against claiming something that cannot be directly seen, but identifiable only 
by its effects, such as a black hole in a binary. 

An interesting recent development is the discovery of "twins" [94]. These au- 
thors show that "Binaries like to be twins"; i.e., something in their birth 
process seems to favor binaries in which the two stars are essentially identical. 
If this is true, it would enhance the double neutron star production, because it 
increases the number of neutron star pairs that are within 4% of each other in 
mass. They did not, however, correct for selection effects, which will probably 
lower somewhat their estimated number of twins. 

Prompted by the paper on twins by Pinsonneault and Stanek [94], Lee et 
al. [66] made an improved calculation of the ratio of neutron star, low-mass 
black hole mergers to double neutron star mergers. The latter could be run 
up if there were lots of twins. However selection effects not included in [94] 
would have to be taken into account before one could estimate the number 
quantitatively. In the non-twin population, Lee et al. found that the number 
of low-mass black-hole, neutron star binaries would be ~ 5 times greater than 
double neutron-star binaries. For detection by LIGO the "chirp" mass would 
be somewhat greater for the former than the latter, because of the greater 
mass for the black holes. 

In summary, it is clear to us that there are ^ 5 times more low-mass black 
hole binaries than double neutron star binaries. As of now, the merging of 
these cannot be distinguished in the short-hard 7-ray bursts. However the 
consequences in terms of the two neutron stars in a binary having to be within 
4% of each other in mass, with small correction made where there has been 
some He red giant mass transfer, are growing in number. We have explained 
why, in the case of SN1987A, the absence of evidence is evidence for a black 
hole. 



4 CONCLUSIONS 

We have been able to arrive at the strangeness condensation phase transition 
critical point ric ^ 3no from three different starting points, the first from the 
zero-density vacuum, the second from the vector manifestation fixed point 
and the third from the Landau Fermi-liquid fixed point. In the two last cases. 
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the estimate was made more reliable by arguments based on renormalization 
group flow which made all of the quantities whose behavior with density was 
unknown rotated out. In particular in the second case, rotated out were the role 
of the scalar condensate (qq), the explicit chiral symmetry breaking through 
T,KN: the role of strange baryons A, E, S, the A(1405) which was often viewed 
as the doorway state for kaon condensation, all of which figure importantly if 
one starts from the T = n = vacuum. We therefore favor the second scenario 
over the others. (In [20] a detailed discussion of these matters is given.) 

We have shown that the hidden local symmetry in the "vector manifesta- 
tion" (HLS-VM) prescribes the medium dependence imposed by (Wilsonian) 
matching to QCD of the quantities entering into the strangeness condensation. 
In particular, since HLS-VM is a fixed point with vanishing gauge coupling 
constant, chiral perturbation theory is well-defined. Assuming that HLS the- 
ory represents QCD near the fixed point, one can efficiently calculate the 
strangeness condensation by fluctuating about the fixed point. It turns out 
that the only degree of freedom we need to take into account is the vector de- 
gree of freedom in the form of on-shell constituent quarks. In other words, the 
strangeness condensation can be calculated in the constituent quark model, 
with the strange quark a passive observer. 

Although we are at a very early stage of understanding, the recent devel- 
opment of holographic QCD, in particular, the Sakai-Sugimoto model, that 
presumably encapsulates QCD at low energy in terms of an infinite tower of 
hidden local fields promises to confirm or invalidate our main theme. Inte- 
grating out the higher members of the tower leaving the lowest could lead 
to HLS-VM theory of Harada and Yamawaki, with the vector manifestation 
representing the effect of higher members of the tower when constrained to 
QCD. At present, it is not known how to go beyond the large and large 
't Hooft limit in the Sakai-Sugimoto model but once one knows how to han- 
dle certain corrections in the dual sector with appropriate matching to gauge 
theory QCD, the characteristic feature of the vector manifestation such as the 
hadronic freedom on which our argument relies could be verified. 

In a neutron star the K~ chemical potential /ik- , essentially the K~ mass, will 
be lowered by the attractive interaction between the quarks in the neutrons 
and protons, chiefly by the former, so that at a certain density, which we 
estimate to be 3no, the fix- becomes equal to fi^-. At that point the electrons 
change into K~ mesons, which go into an S-wave Bose condensate. With the 
accompanying lowering of pressure, the star drops into a black hole. This 
we believe is the flrst phase transition dense compact star matter undergoes 
which determines the fate of star, that is whether it becomes a neutron star 
or a black hole. 

We believe that this scenario describes the fate of the supernova explosion 
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SN1987A. What we are sTirc of is that in the evolution of binary neutron 
stars, the first-born neutron star will go into a black hole if the two progenitor 
giants are not within 4% of each other in mass. (For the low-mass neutron-star 
binary the situation is slightly shifted, as we explained in the text.) 

Since our main early interest in kaon condensation was in the disappearance 
of evidences of the existence of neutron star in SN1987A after it had emitted 
neutrinos for about 12 seconds, we felt it fitting to adduce even more firm 
astrophysical observations that shows that there are a large number of low- 
mass black holes in the Universe. Most specific is that we showed that double 
neutron stars in binaries must be within 4% of each other in mass; otherwise 
the binary would end up being one of a neutron star and a low-mass black 
hole. Observations support the necessity of the neutron star being within 4% 
of each other in mass within any given binary. This also means that there 
are ~ 5 times more neutron-star, low-mass black hole binaries than double 
neutron-star binaries in the Galaxy. 

Finally, putting together all of the above leads us back to the proposal of 
Brown and Bethe [37] that the maximum neutron star mass is ~ 1.5Mq, not 
much higher than the mass of the Hulse- Taylor pulsar 1913-1-17 of 1.44Mq. 
This is the only accurately measured mass of the more massive neutron stars. 
We believe, for example, the probability that J0751 + 1807, of mass 2.1±0.2 but 
with 95% probability > l.GM© has mass ~ l.SM© to be substantially higher 
than for the spectrum of neutron stars in binaries that would exist were 2.IM0 
neutron stars to be stable. The announced reduction of the measured mass 
of J0751-I-1807 to ~ I.26M0 goes in the direction consistent with the above 
reasoning. 

The Harada and Yamawaki HLS-VM had predictions that countered well es- 
tablished behefs. Perhaps the most controversial was the prediction in HLS- 
VM of hadronic freedom, that the interactions between hadrons go to zero 
as the temperature goes up to from below and the density goes up to 
from no. This directly contradicted the accepted belief in heavy- ion commu- 
nity that interactions between hadrons became more plentiful and stronger 
as T goes up to from below. We managed to find what we beheve to be 
the crucial experiment to decide this issue, namely the STAR experiment of 
peripheral Au+Au collisions which measured ratio [26]. We discussed 

this experiment briefly in Sec. 5. The experiment was so useful because in 
this peripheral experiment the freeze out temperature was equal to the flash 
temperature, the temperature at which the hadrons went on shell. The results 
of this STAR experiment directly supported the concept of hadronic freedom 
in temperature. We are proposing that the same is the case in density and 
applies to kaon condensation. We found in this paper that HLS-VM gave un- 
ambiguous statements about how to handle the medium dependencies and 
that is why we are so sure of the results. Needless to say, the crucial element 
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in this argument is the existence of the vector manifestation (VM) fixed point 
which requires the vanishing of the vector meson masses (in the chiral hmit). 
The VM in Harada-Yamawaki HLS theory is an unambiguously falsifiable pre- 
diction: it could be invalidated if an unquenched lattice measurement of the 
spectral function of the p channel unquestionably showed that the (real-time) 
mass of the p meson does not vanish in the chiral limit at T = T^. 

The question we have not addressed here is this: if kaon condensation is in- 
evitable in compact stars once density reaches ~ Stiq, will other forms of exotic 
states that can take place at higher densities, say, > ~ SriQ, such as color 
superconductivity, color flavor locking etc. be relevant for the physics of com- 
pact stars? We can offer no definitive answer to this question. Our conjecture 
is that they will be hidden in the "black-hole nothingness" inaccessible to 
the observers. The question can perhaps be answered when holographic QCD 
becomes sufficiently sophisticated and predictive. 
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